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Abstract 

We present a study of a water-in-oil microemulsion in which surfactant coated water nanodroplets 
are dispersed in the isotropic phase of the thermotropic liquid crystal 5CB. As the temperature is 
lowered below the isotropic to nematic phase transition of pure 5CB, the system displays a demixing 
transition leading to a coexistence of a droplet rich isotropic phase with a droplet poor nematic. 
The transition is anticipated, in the high T side, by increasing pretransitional fluctuations in 5CB 
molecular orientation and in the nanodroplet concentration. The observed phase behavior supports 
the notion that the nanosized droplets, while large enough for their statistical behavior to be probed 
via light scattering, are also small enough to act as impurities, disturbing the local orientational 
ordering of the liquid crystal and thus experiencing pretransitional attractive interaction mediated 
by paranematic fluctuations. The pretransitional behavior, together with the topology of the phase 
diagram, can be understood on the basis of a diluted Lebwohl-Lasher model which describes the 
nanodroplets simply as holes in the liquid crystal. 
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I. INTRODUCTION 



Liquid Crystals (LC), because of their soft spontaneous ordering, have been found to 
produce interesting properties when used in the production of micro- or nano-structured 



heterogenous materials 



l|. Relevant examples are given 



by the large electro-optic suscep- 
2]), by the novel thermodynamic 



tibility of LC confined in a solid matrix (e.g. PDLC 
properties of LC in random solid networks (e.g. LC in silica aerogels jj, |J|), by the un- 
precedented interparticle forces observed between LC suspended micron-sized colloids (e.g. 
surfactant stabilized water emulsions 5|). In these three geometries (confinement, bicontin- 
uous structures, LC as suspending medium), the organization is surface dominated. Given 
the soft LC ordering, forces arising at the interfaces effectively compete with bulk forces, 
yielding a wide combination of simultaneous effects such as wetting and pre-wetting, critical 
fluctuations, topological defects, surface phase transitions and elastic deformation. 

Within this framework, a new intriguing possibility has been recently proposed: materials 
combining the tendency to molecular ordering of thermotropic LC, with the tendency for 
self aggregation of amphiphilic molecules in solution. These systems are simply obtained by 
replacing the oil phase of surfactant-water-oil mixtures with a thermotropic LC. When the 
molecular structures of the constituents are chosen in an appropriate manner, these mixtures 
organize in ways still largely unexplored. 



After the pioneering work by Marcus 



6|, water-surfactant-LC mixtures have been exten- 



sively studied first by Yamamoto and Tanaka 



7( , and later by us 



Yamamoto and Tanaka 



reported the "transparent nematic" phase in a composite nematic liquid crystal (penthyl- 
cyanobiphenyl, 5CB), double tailed ionic surfactant (didodecyl-dimethyl-ammonium bro- 
mide, DDAB), and water system. The surfactant, which tends to form interfaces concave 
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toward the water at organic-water interfaces, forms a water-in-oil type microemulsion of 
inverse micelles in the 5CB host. Upon cooling, these mixtures were reported to exhibit 
two phase transitions as indicated by double peaked DSC scans: first to a new optically 
isotropic transparent nematic phase and then to the coexistence of bulklike nematic (N) 
and surfactant-rich isotropic (I) phases. According to the interpretation, the new phase is 
actually a nematic phase in which elasticity and radial boundary conditions on the micelle 
surface combine to produce a nematic structure where the director is distorted on the lOnm 
length scale. This picture appears quite similar to that offered by the behavior of micron- 
sized spheres in nematic media, implying that micellar nano-inclusions, despite their size, 
comparable to the LC molecules, can also be accounted for using concepts typically suited for 
larger length scales, such as use of the molecular director, continuous elastic distortions and 
surface anchoring. However, the additional results provided by our investigation on the same 



system 



8| contradict the notion of a "transparent nematic" phase and support a different 



interpretation: upon cooling from an isotropic phase resembling an ordinary microemulsion 
of inverse micelles, the system undergoes a demixing transition in which the micelle-poor 
phase develops long range nematic order. This phase behavior, already observed 9| and 
predicted 10| for molecular scale contaminants such as polymers, can be understood on 
the basis of a very simple model, in which micelles are considered as "holes" in the LC. 
Experiments and theory thus combine to picture this system - LC with nano-inclusions - 
in a way radically different from its analog on the micron scale: nano-scale impurities act 
much more as molecular contaminants - lowering the local degree of LC ordering - than as 
microparticles - characterized by specific surface anchoring of the surrounding LC. More- 
over we observe, in agreement with theoretical predictions, the presence of pretransitional 
micellar density fluctuations anticipating, on the high temperature side, the demixing. Such 



fluctuations can be understood as resulting from intermicellar attraction, in turn depending 
on the orientational paranematic fluctuation of the LC matrix. Paranematic fluctuation me- 
diated interactions have never been observed before, although the presence of Casimir-style 
interactions is generally expected. However, because of the first order nature of the I-N 
transition of thermotropic LC, the correlation lengths developed in the I phase extend only 
up to 10 — 20nm, making pseudo-Casimir forces quite difficult to observe. 

In this article we present a new body of results on the same microemulsion (hereafter 
referred to as 5DH jiem) described in Ref. js|. New data obtained by Electric Birefringence 
Spectroscopy (EBS) and by static and dynamic Light Scattering (LS), ( Sections IIII Bl and 
IIII C|) coherently match the description of the system previously offered 8[, as explained 
in the discussion (Section^). In Section HV| we present a theoretical analysis of a simple 
spin model, which has been extended in comparison with the one described in Ref. 8| 
in few respects, as further elaborated in Section IIV1 The relevance of the model in the 
interpretation of the experiment is discussed in Section In Section IS we also compare the 
behavior of 5DHfiem - in which the inverse micelles can be viewed as providing annealed 
disorder to the LC - with the behavior of LC incorporating quenched disorder as described 
in a wide body of previous work. 



II. EXPERIMENTS 

A. Sample preparation 

The system studied in this paper is a three component mixture, whose components 
are 5CB (penthyl-cyanobiphenyl, M w = 249.36) from Merck, DDAB (didodecyl-dimethyl- 
ammonium bromide, M w = 462.65) from Sigma-Aldrich and de-ionized distilled water. Sam- 
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pies were prepared by adding 5CB to a well homogenized mixture of DDAB and water at 
room temperature. The samples were then hermetically closed, kept at ~ 60°C and stirred 
to ensure homogeneity. 

The system behavior, and particularly the transition temperature, was found to be very 
sensitive to the concentration of water and thus reliable sample preparation required careful 
control of evaporation and hydration of components and containers. No other critical matter 
resulted emerged in the preparation procedure. At 60°C, the system was indeed found to 
easily mix and promptly homogenize, yielding a uniform transparent isotropic state. 

All the data presented in this article will refer to sample on the dilution line characterized 
by weight ratio viiddabI^h^o — 9. We will indicate with the weight ratio <fi = (, mr J^° mo ) , 
where rri£, = itiddab + rnn 2 o is the mass of the LC dispersed material. Being pddab ~ 
S70kg/m 3 and p§cB ~ 1022kg/m 3 we can calculate the volume fraction of dispersed material 
4>v ~ 4> ' 1-17, and the overall density of the dispersed material pr> ~ 881kg/m 3 . 

As it will be shown below, dispersed water and surfactant aggregate into inverted micelles 
having a water core, whose number density N m i C can be calculated once the volume of the 
single droplet (v m i C ) is known. 

We prepared different set of samples for the different experimental investigations i.e. de- 
termination of phase diagram, Kerr constant measurements and light scattering experiments. 
We have found that different samples prepared with the same method display phase separa- 
tion at almost identical temperatures, which confirms that the adopted sample preparation 
procedure is reliable. 
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B. Volumetric determination of the phase diagram 

For temperature larger than the bulk T^j, all the produced 5DHfiem samples are opti- 
cally isotropic. As T is lowered, a demixing phase transition takes place at a dependent 
temperature Tdm(<P)- At T < Tdm{4>) samples separate into an isotropic phase coexisting 
with a denser optically uniaxial (nematic) phase. By maintaining the samples at constant 
temperature and waiting several hours, the two coexisting phases macroscopically separate 
under the effect of gravity. The samples appear as shown in Fig. Q ( a ); where the highly 
scattering bottom nematic phase appears white. The two phases are easily distinguish- 
able and the relative volume can be measured by taking pictures and analyzing the digital 
images. In order to extract from volumetric measurements the phase boundary of the coex- 
istence region, we prepared - by successive dilution - four different samples: s\ (0 = 0.02), 
s 2 {4> = 0.05), s 3 (<ft = 0.075), s 4 (<f) = 0.15). These samples were kept into a thermostatic 
chamber with a stability of ±5 • 10~ 3 i^ and observed by a digital camera at various temper- 
atures. In the inset of Fig. |2]we show the T dependence of the ratio, Vi/Vtot, between the 
volume of the isotropic phase and the total volume of the sample. Sample concentrations 
were chosen so that, in a wide temperature interval, at least two samples simultaneously 
have separated phases large enough to be easily measured (typically larger than 10% of the 
total volume). At each temperature we have extracted the micellar weight ratios </> 7 (T) and 
4> N (T), enclosing the coexistence region of the 5DHfiem, shown as dots in Fig. |21 

In this analysis we have assumed that the behavior of 5DHfiem can be understood as a 
pseudo-single component system, i.e. that micellar aggregates don't change their structure 
as the phase transition takes place, but simply redistribute between the two phases. This 
assumption is confirmed by light scattering measurements in the isotropic phase, as discussed 
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in the next Sections. Accordingly, from the standpoint of micelles, the phase transition is 
a simple demixing in which they preferentially accumulate in the isotropic phase. Since 
micelles are less dense than 5CB, their presence increases the density difference between 
isotropic and nematic phase. 

C. Electric Birefringence Spectroscopy 

Measurement of the Kerr coefficient have been performed by applying pulses of ac voltage. 
The pulse duration (0.1 — 10ms) was abundantly enough for the system to reach the steady 
state. The frequency v of the applied field ranges from 1 kHz to 300 MHz, the upper limit 
being set by inductive loops introducing distortions in the detection of the field applied to the 
cell electrodes. For this reason the wiring has been designed to keep wire loops to a minimum. 
The field induced birefringence is, at steady state, composed by a dc component and by an 
ac component of frequency 2u. While the former is related solely to the amplitude of the 
local dielectric anisotropy characterizing the nematic fluctuations in the isotropic phase, the 
latter is instead strongly dominated by the dynamics of the paranematic fluctuations. The 
set up used for these measurements is not equipped with a fast enough detection optics 
and transient digitizer to reliably measure the amplitude of the ac component. It is instead 
very accurate in the measurement of the dc component: the use of a quarter- wave plate in 
between cell and analyzer, the automated set-up enabling measurements at various angular 
positions of the analyzer and various field amplitudes, as well as the accumulation of a large 
number of pulse response at each studied frequency, allows extracting the field induced non 
oscillatory component of the birefringence An^c = nn — n±_ with an error of few percent. 
The measured An dc , has been found, as expected, to be proportional to the square of the 
electric field amplitude E, hence enabling a meaningful extraction of the Kerr constant B{y) 
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defined as: 



where A = 633nm is the optical wavelength used in the experiment. In Figs. El and 0] we 
show B(y) spectra measured in bulk 5CB and in a 5DHfiem sample. As clearly visible, our 
experimental frequency range includes the main features of the spectra. 

D. Static and Dynamic Light Scattering 

Light scattering measurements have been performed by a highly automated experimental 
set up, featuring a He-Ne laser source, a temperature control with a 5 • 10~ 3 K stability 
and collection optics coupling scattered light into optical fibers. We used multimode fibers 
to measure the intensity of both polarized and depolarized scattered light (Ivv and Ivh 
respectively). Single mode fibers were instead used to measure intensity autocorrelation 
functions of the polarized and depolarized scattered light (G 2 yv and G 2 yH, respectively). 
Use of fiber beam-splitter and cross correlation has enabled to reliably extract correlation 
functions down to a retardation time of 25 ns, a limit set by the electronics of the BI9000 cor- 
relator board used in the experiment. All static and dynamic measurements were performed 
at fixed angle 9 = 90°. 

In Figs. I5|fil we show the polarized and depolarized scattered intensity for bulk 5CB and 
5DH[iem samples. While Ivv and Ivh measured in each sample at various T are on the 
same intensity scale, comparison between intensities scattered by different samples should 
be considered only as semi-quantitative since cell replacement required realignment of the 
collection optics. 
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Fig. [7| shows typical Giyv an d G^vh measured by studying bulk 5DHfiem samples. 
The inset of Fig. [7| instead compares G2,vv as measured in = 0.075 and in <fi = 0.15 
5DHfiem. 



III. DESCRIPTION AND ANALYSIS OF THE RESULTS 
A. The phase diagram 

The most relevant features of the phase diagram in Fig. El are the coexistence of a 
micelle poor nematic and a micelle rich isotropic over an extended temperature range and 
the presence of a re-entrant nematic phase. 

It is known that inclusion of molecular dopants in a nematogenic compound results in 
extending the isotropic-nematic phase coexistence - confined to T = Tni in bulk LC - to 
a wide T interval. This phenomenon has been experimentally studied in detail for low 



ll| . The authors also 



contaminant dilution of aliphatic molecules having various shape 
present a thermodynamic analysis of the T^j((p) slopes of the isotropic and nematic phase 
boundaries (approximated by straight lines converging to T/vj). In particular they showed 
that nematic-isotropic thermodynamic equilibrium implies that 



-1 - 1 = (2) 

(3 N (3 1 R 

dT R 

where (3^,1 = are the are the slopes of the (respectively nematic and isotropic) phase 
boundaries enclosing the coexistence region, T R = x is the contaminant concentration 
expressed in mole fraction, R is the gas constant and ASni is the N-I transition entropy for 
Dulk 5CB that can be extracted from calorimetric measurements, namely AS^j/R ~ 0.25 
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According to Eq. E the quantity (fi^ 1 — (3j l ) is independent from the nature of the 
contaminant. Ref. ll| reports that, for molecules of elongated shape, (3^ ~ —0.58 and 
f3j ~ —0.51, thus approximately confirming the predicted behavior. A successive study of 
the effect on the I-N transition of 5CB dope d by surfactant molecules of the polyoxethylene 



family reports /3n ~ —1.35 and (3x ~ —0.95 



121 ] . also is agreement with Eq. El 



The extrapolated initial slopes of the phase diagram in Fig. Hare: dTft/d<f> 1.23 and 

dT^/dcj) ~ —0.20. By assuming that all contaminants have the molecular weight of DDAB 
(i.e. disregarding the presence of water), we obtain /3n ~ —2.6 and (3j ~ —0.35, while, if 
we calculate exactly the total molar fraction of contaminants (water included) we obtain 
Pn ~ —0.75 and [3i ~ —0.10. In both cases the prediction of Eq. |2]is not confirmed since 
(Pn 1 ~ P7 1 ) i s 10 to 30 times larger than predicted. Indeed, the bDHfiem phase diagram 
is evidently characterized by a separation between the isotropic and the nematic phase 



ll| and 



12j. This analysis suggests that, in 



boundaries much larger than those found in refs 
calculating j3, we are overestimating the molar concentration of dopants. Indeed, completely 
different estimates are obtained if one considers, as contaminants, not the single monomeric 
dispersed surfactant molecules, but the micellar aggregates, as discussed in section |V] after 
micellar sizes and interactions are introduced. 

Although predicted by many previous theoretical analysis, the presence of a re-entrant 
nematic phase has been experimentally observed only quite recently 
of the liquid crystal E7 and silicone oil. The nematic phase has been found to extend to a 
maximum contaminant concentration of 3 % vol. at about 10°C below TJvj, i.e. similar but 
slightly larger than the one in Fig. [21 The physical origin of reentrant nematic phase can 
be sought in the combination of the positive T vs. nematic phase boundary expected at 
the lowest temperatures with the negative T vs. <fi phase boundary resulting predicted for 
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9( by studying mixtures 



temperatures close to the NI transition. While at the limit of low temperature the nematic 
crystal expels all impurities, as T is increased the entropic gain for impurity solubilization 
determines an increment of the maximum micellar concentration in the N phase. Upon 
increasing further the temperature, the free energy difference between the bulk N and I 
phase is reduced, vanishing at Tjv/- Therefore for each density of impurities solubilized 
in the nematic at low T, a maximum T is found where the free energy of the nematic 
with impurities equals the energy of the isotropic phase. At this temperature it pays for 
the system to phase separate into two phases as described by equation |21 The cross-over 
between the entropically driven impurity solubilization at low T and the impurity driven 
phase instability at high T leads to the observed reentrant behavior. 

Phase diagrams similar to those in Fig. |2]have also been found in a somewhat different 
research contexts. The contaminant induced I-N phase separation is exploited to produce 
Polymer Dispersed Liquid Crystal materials. An isotropic mixture of isotropic LC and 
polymers is destabilized by a T quench. The isotropic phase hosting nematic nuclei is 
subsequently hardened by photo or chemically induced cross polymer linking. Such a phase 
separation has therefore received particular attention. Its experimental characterization, 



performed by studying a mixture of E7 and polybenzylmethacrylate [10 1, is qualitatively 
similar to the isotropic phase boundary in Fig. |2] although few experimental points may 
suggest the presence of an additional isotropic-isotropic phase separation. 

B. Electric birefringence results 

The onset of local nematic structuring, as the one predicted by Yamamoto and Tanaka's 
"transparent nematic", would quite likely result in an electro-optic susceptibility larger than 
the conventional isotropic phase, as also suggested by the reported shear birefringence of such 
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a phase. This is why measurements of electro-optic response by EBS appear particularly 
appropriate to study the 5DWfj,em. 



1. Bulk 5CB 



According to the mean field description, the Kerr coefficient for bulk LC is given by 



13] 



_ e Ae*(z/)An 
B[V) ~ 3Xa(T - T*j (3) 

where a, Ae^ and An are, respectively, the first coefficient in the Landau-DeGennes 
expansion and the dielectric anisotropy and birefringence of a completely ordered nematic. 
Ae* = ej^ — e^, the parallel and perpendicular component of the dielectric coefficient of a 
perfectly ordered nematic (the superscript X refer to nematic crystal), gauges the ampli- 
tude of the coupling of E with the induced nematic order, and is the only term in Eq. |3] 
depending on the frequency of the applied electric field. However, the contribution of to 
the frequency dependence of B is negligible: the permanent dipole of 5CB is aligned with 
the molecule and thus only electronic polarization contributes to e^, which is therefore not 
expected to have observable frequency dependence in our experimental EBS spectra. The 
dominant contribution of arises from molecular reorientation, its characteristic frequency 
deriving from the kinetics of molecular reorientation around the short axis. We thus expect 
EBS spectra for isotropic bulk 5CB to be close to its dielectric spectra. 

In short, Eq. El indicates that EBS spectra enable accessing two independent quantities: 
the dynamics of the individual molecules from the frequency dependence of B(u), and the 
degree of local paranematic ordering from the low frequency value of B. It is worth noticing 
that this last quantity cannot be detected by ordinary dielectric spectroscopy. 
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With respect to the more conventional pulsed electric birefringence, EBS better enables 
a distinction between the various processes taking place in the different frequency ranges. 
This is useful when low frequency artifacts, such as electrode polarization and space charge 
accumulations, typically frequency dependent, are present. In practice, our data, shown in 
Fig. El show little frequency dependence in the kHz regime. We have fitted, for v > 1MHz 
the measured B(y) with the sum of a frequency independent value, B e , expressing the Kerr 
coefficient due to the anisotropy in the electronic polarization of 5CB molecules, and of a 
Debye type relaxation describing the frequency dependence of the polarization due to the 
reorientation of permanent dipoles: 

BM = ft + IT(l^ (4) 

In the fitting procedure, B e , Bd and v c are used as fitting parameters. The good quality 
of the fits indicates that a single relaxation frequency is enough to describe the whole polar- 
ization process. This may appear surprising since previous investigations by high frequency 
and time-domain dielectric spectroscopy on cianobiphenyls in the isotropic phase showed 



evident deviations from Debye behavior 1J| 1£ 



lfi| . In fact, Debye behavior is instead 



regularly observed when studying the ey spectra in the nematic phase. We argue that, since 
B{y) reflects the frequency dependence of ejF, it is not actually unexpected that it behaves 
more similarly to the nematic en than to the isotropic e. This phenomenology together with 



more data on cianobiphenyls compounds, will be discussed in a forthcoming publication XL \ ■ 
The three fitting parameters are found to depend on temperature. At all temperatures, 
we find that B e ~ B^/30. This ratio is in rough agreement with what expected by comparing 



dielectric (y < u c ) and optical (v > u c ) anisotropies: from literature values 



2 1, in the nematic 



phase (njj — n 2 ± )/(e\\ — e±) ~ 1/20. The factor 1.5 in between estimate is probably due to 
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the fact that our experiments are performed in the I phase, at values of field induce nematic 
order parameter much smaller than achievable in the N phase (S(E = lOV/mm) ~ 1 • 10~ 6 ) 
where the above ratio has instead been evaluated. The temperature dependence of B e and 
Bd is shown in Fig. |HJ As indicated, the uncertainty on B e is much larger than for B^, and 
its T dependence has to be considered only qualitatively. In the same figure we have plotted 
the fit of B d (T) to A 5CB (T-T*y\ yielding A 5CB = 1.91(r 1( W- 2 and T* = 33.4±0.1°C. 
Both the quality of the fit and the value of T* in agreement to previous determinations, 
confirms that the amplitude of the Kerr effect is well described by a simple mean field 
theory. 

In Fig. El we plot v c as a function of 1/T. The Arrhenius-like behavior, 



v c oc e (" A ™ (5) 

is evident, where AH is the activation energy. We find AH = 41.0 kJ/mol. This 
number should be compared with AH = 33 kJ/mol, obtained by fitting the Arrhenius 



dependence of the cut-off frequency of the 5CB dielectric spectra in the isotropic phase 



id- 



According to mean-field expression, both EBS and dielectric characteristic frequencies should 



be proportional to the fluid viscosity 



191 ] . Such proportionality is partially confirmed when 



AH t is compared with activation energy extracted from by the Arrhenius T dependence of 



5CB viscosity: from literature we get: AH V = 29.0 kJ/mol 



2l|, AH V = 28.3 kJ/mol 



, AH V = 44.9 kJ/mol 



22j. Despite the rather surprising numerical discrepancies, 
it is very clear from the data that EBS, in analogy to viscosity, probes the dynamical 
response of individual molecules, undergoing activated processes of angular reorientation. 
We will exploit the access to local molecular dynamics through EBS to evaluate how much 
the micellar inclusions affect the local viscosity, a notion necessary to interpret micellar 
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diffusion. 



2. bDHnem 

As a simple inspection of the data in Fig. |H] immediately reveals, EBS results for the 
DDAB, water and 5CB mixtures are qualitatively similar to the bulk data although depressed 
in temperature and limited in amplitude. Although this is already enough to make clear that 
no extravagant molecular organization is present in the system, it is nevertheless interesting 
to analyze the results as it has been done for bulk 5CB. 

All EBS spectra can be fitted with Debye dispersions (Eq. HJ as in the bulk. In Fig. Owe 
show the low frequency Kerr coefficient Bd as a function of temperature. The lowest value 
achieved for each given (0 = 0.15 and = 0.075) corresponds to the temperature at which 
demixing occurs, T DM ((p = 0.075) ~ 31°C and T DM ((p = 0.15) ~ 27.5°C In agreement with 



visual observations ( section |IIB|) . the isotropic phase of 5DH/j,em extends to temperatures 
lower than the bulk T/vj. By fitting Bd to A^cb{^- — <P)/{T — T*), we obtain a value for 
T* depending on 0, namely T*(0 = 0.075) = 27.3°C and T*(0 = 0.15) = 21.8°C. As 
increases, the maximum value attained by the Kerr constant decreases, in agreement with 
the increased difference between Tdm and T*. These data, which will be discussed in Section 
1111 C 2l together with the light scattering results, indicate that the electro-optic susceptibility 
of the mixture is dominated by the LC component and that its behavior is simply the one 
expected for an isotropic LC approaching the nematic transition with shifted T*. 

Since the individual molecular dynamics is typically unaffected by the onset of parane- 
matic fluctuations 23|, as also confirmed by our data for bulk 5CB, it is expected that the 
characteristic frequencies extracted from EBS 5DHfiem would also not be influenced by 
pre-demixing phenomena. In Fig. Elwe plot v c vs. 1 /T for the comparison with bulk results. 
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Within the experimental uncertainty shown by the spread in the data point, the Arrhenius 
behavior of the 5DHfiem v c data coincides with the bulk one, indicating that the activation 
energy for molecular reorientation is equal in the two systems. Thus the presence of the 
micelles doesn't affect the molecular environment experienced by the LC molecules, which 
remains basically isotropic-like down to Tdm- Given the large fraction of 5CB molecules in 
contact with the DDAB molecules around the micelles (an estimate will be given in Section 
EJ), this result clearly indicates that no significant anchoring of LC molecules is induced at 
the micelle-LC interface. 

In Fig. EI] we show the characteristic time r c = (27ru c )~ 1 extracted from the data in Fig. 
HJ to be compared with the characteristic time for paranematic fluctuations, as discussed 
below. 



C. Light Scattering Results: Extracting Paranematic and Micellar Contributions 

1. Bulk 5CB 

Light scattering provides valuable information on both the statics and dynamics of pre- 
transitional paranematic fluctuations through, respectively, the intensity and the correlation 
function of the scattered light. Because of the first order nature of the I-N phase transi- 
tion, light scattering results are usually accurately described by the mean field theories. 
Accordingly, the scattered intensity, as described by the Rayleigh ratio R, grows as 



R ~ a(T- T*) (6) 
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while the correlation time of the paranematic fluctuations is given by 



23] 



_ k B Tr] 

TPN ~ a(T - T*) (7) 

where rj is the viscosity. 

Our bulk data are in agreement with the mean field description and with previous inves- 
tigations. Measured lyy and Ivh intensities are well fit by EqEJwith added a background 
constant, accounting for stray light. Fig. |B1 shows such a fit for Ivh, yielding T* = 33.8°C 
(and thus Tni — T* = 0.9°C), in acceptable agreement with literature [24| and EBS data 



analysis. Polarized and depolarized intensity are found to be proportional, as predicted 23] : 



{Ivh) = \{Ivv) (8) 
This is shown in Fig. |3] where cross symbols represent | • (Ivh) an d nicely overlay on 
(Ivv) data. 

The correlation functions obtained in the bulk can be well fit with exponential decays. 
The analysis of the VV correlation functions yields to values of tpn (see Fig. EH full dots) 
about 15 % larger than those obtained from VH scattering. This observation is in line with 
prediction in the frame of mean field analysis as a consequence of coupling between shear 
viscosity and density fluctuations, which bears detectable effects in the VH Rayleigh spectra 
. On the basis of equation (23) of Ref. 



and by using values obtained from experiments 
on MBBA, we would predict tpn,vv/tpn,vh ~ 1-3, about double of what observed in 5CB. 

t pn/v can be well approximated by (T — T*)" 1 , as predicted by Eq. [7| and shown in 
Fig. m The extracted divergence paranematic temperature is T* = 34.4°C is too high with 
respect to the one obtained from Ipn in the same light scattering experiment. The origin of 



this difference appears to be seeded in the behavior of r\ 
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2. 5DH/iem 

Microemulsion data require some additional analysis. The measured correlation functions 
(FigEJ) very clearly demonstrate the presence of two distinct processes: a fast one, appearing 
in both VV and VH correlations, clearly due to paranematic fluctuations, and a slower 
one, appearing only in the VV correlations, detectable only in the microemulsions. This 
last component, indicating the presence of aggregates within the 5CB matrix, is easily 
interpreted as intensity scattered by the inverted micelles. Thus the total scattered intensity 
Itot — Ipn + Im is the incoherent sum of intensity scattered by LC orientational fluctuation 
(Ipn) and by micelles (Im)- Specifically 

Iyy = ItoT,VV = IpN,VV + hi (9) 

while 

IvH = ItOT,VH = IpN,VH (10) 

The large difference between the correlation times of the two components makes it possible 
to identify them. Let's indicate with T m ^ ~ 1 — 10/xs a correlation time long with respect to 
paranematic fluctuation time and short with respect to micellar diffusional times. Because 
of the statistical independence of the two terms we have 

G2,vv(^) ~ G 2 yv(Tmid) = c ((IpN,Vv) + 2(IpN,w) (hi)) 
G2,w(T m id) - G 2 yv(oo) = c{I 2 M ) 

where c is a factor related to the collection optics. Combining the two equations, (Im) can 
be derived from the total polarized scattered intensity and from the decay of the correlation 
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functions. Alternatively, (Im) can be extracted from the polarized scattered intensity (Ivv) 
as 

(Im) = (Ivv) - Ivh) (12) 

where we have used the notion that (IpNyv) — 4/3(Ivh) as observed for paranematic 
scattering in bulk 5CB (eq. |BJ). Fig. ^Jcompares the values of (Im) vs. T obtained from the 
two procedures. The comparison clearly confirms that the T dependence of (Im) is indeed 
reliably extracted. The determination of (Im) from DLS data, however, requires extracting 
Cr2,w(0) from the correlation functions. This has been done by fitting Gzyvi^) by a sum of 
two exponentials (see Fig. EJ) enabling determining G ! 2,w(0), as well as tpn and the micellar 
correlation time tm- tpn obtained from G2,vh(t) are more reliable than those obtained from 
G 2 yv(T) since in this last case the intensity scattered by the micelles, basically constant on 
the times scale of paranematic fluctuations, beats with the light diffused by orientational 
fluctuations thus introducing an heterodyne contribution to the correlation function. There- 
fore, tpn extracted from G2,vv(j) would, in principle, require a compensating correction, 
unnecessary in the VH scattering. Although in practice such a correction is small since 
Im/Ilc < 15% always, we will in the following refer to r PN obtained from G2.yH(j). 

In summary, light scattering measurements enable extracting (Ipn) , (Im) , tpn, tm as a 
function of T and for the different samples. In what follows we will separately discuss results 
pertaining the two sides of the 5DHnem behavior: LC in the presence of nanoinclusions 
and micelles within a solvent with ordering fluctuations. 

Inspection of (Ipn) vs. T measured in the 5DHfiem shown in Fig. |U1 reveals that: (i) at 
any given T, paranematic fluctuations are depressed by the presence of the micelles; (ii) an 
isotropic phase like behavior extends down to Tp>Mi i- e - even at T lower than the bulk Tjvj; 
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(iii) the maximum value of (Ipn) for each given preparation, obtained for T = Tdm, is a 
decreasing function of </>. 

By fitting I PN (T) as for the bulk data we obtain T*(0 = 0.075) = 27.2°C and T*(0 = 
0.15) = 21.9°C, closely matching the results obtained from EBS. These results are reported 
in the phase diagram in Fig. El where the dotted line interpolating the data guesses the 
behavior of T* (</>). 

tpn/v data vs. T are shown in Fig. ^2 Given the uncertainty on the T dependence of 
the viscosity, as discussed above, we prefer not to fit the data, but to evaluate the magnitude 
of the temperature shift induced by the presence of the micelles by finding the best overlap 
between the 5DHfiem data and the bulk 5CB data shifted down in temperature of 6T as 
shown in Fig. ^2 Best overlaps have been obtained with 5T = 6.25°C(0 = 0.075) and 5T = 
ll.O o C(0 = 0.15), to be compared, respectively, with T* (bulk) -T* ((f) = 0.075) = 6.6°C and 
with T*(bulk) - T*(<f> = 0.075) = 11.9°C as obtained in static LS. 

The overall agreement between the T* from static and dynamic LS and from the low 
frequency Kerr measurements, supports the notion that the annealed disordering provided 
by the micelles is affecting the LC behavior simply by shifting it to lower temperatures and 
by increasingly precluding the development of paranematic fluctuations by the appearance 
of the demixing transition, i.e. T DM ((j) = 0.15) - T*(</> = 0.15) [~ 10°C] > Tdm (0 = 
0.075) - T*(0 = 0.075)[~ 6°C] > T NI - T*(bulk)[~ 1°C\. 

3. Micellar behavior 

The temperature dependence of (Im) for the two concentrations is shown in Fig. The 
most interesting characteristic of these data is the growth of (Im) as T decreases, which 
cannot be understood simply in terms of T dependence of refractive indices, which, in the 
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T interval of Fig. vary at most by 1 %. Analogous deviation from simple behavior 
is observed in the diffusion dynamics of the micelles. The micellar diffusion coefficient 
Dm extracted from r M x = 2DmQ 2 (where q is the scattering vector) is shown in Fig. EH 
and compared with the Dm which would be expected from the Einstein equation Dm = 
(k B T) I '(6irrjR) in the case of spheres of radius R = 18A diffusing in the isotropic phase, 
where the viscosity is taken from Ref. j2l|. The radius has been chosen so to approximate 
the observed Dm at the highest measured temperatures. As evident from the comparison, in 
the proximity of the demixing transition, the measured dynamics is slower than the expected 
one by a factor of about 2 in both the = 0.15 and the = 0.075 samples. 

The comparison between Dm{T,4>) and (ksT) /(6irRr](T)) is based on the assumption 
that the local viscosity experienced by inverted micelles in the 5DHfiem is unchanged with 
respect to bulk 5CB and that, for T < T/vj, where no bulk data for isotropic viscosity exists, 
T}(T) is given by the same Arrhenius dependence as in the bulk. This assumption is is based 
on our EBS results on the Kerr z/ c , which proved that the local molecular environment in the 
isotropic phase of bDH^iem overlaps the low T extrapolation of the bulk Arrhenius. Given 



that viscosity and molecular diffusion share their basic physical origin |25(, we conclude 
that similar T dependence is expected for r)(T). Use of another among the published 5CB 
viscosities would not change the comparison in Fig. ED Namely, by using r/(T) as in Ref. 
we would obtain R = 21A and a maximum ratio between expected and measured D of 



about 2.1, while by using rj(T) as in Ref. 



we would obtain R = 24 A and a maximum 



ratio between expected and measured Dm of about 1.8. Given this dispersion of the R 
values, in the following discussion we will assume R = 20 ± 2A. 

According to our data, and within our experimental error, the diffusion coefficients mea- 
sured for the two 5DHfiem micellar concentrations vs. T are basically the same (see Fig. 
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[T3*j) . Although this result may appear as a straightforward indication of having droplets of 
the same size in the two samples, some analysis is instead required. In fact, = 0.15 is a 
rather significant volume fraction, at which effect from interparticle interactions may be ex- 



pected. When expressed through the virial expansion 
measured at a given particle concentration is 



the collective diffusion coefficient 



D M = D Mfi {l + k D <f>) (13) 

where Dm,o is the diffusion coefficient at infinite dilution and kn is the virial coefficient 
expressing the correction to collective diffusion arising from interparticle interactions. Ac- 
cording to the data here presented, the product ko(fi assumes the same value in the two 
different samples, thus indicating a independent D, and possibly kr, = 0. 

In order to check the robustness of this evidence and also to explore possible deviations 
from simple one-component Maxwell rule for the phase coexistence, we have studied the 
upper phase - denser in micelles - obtained by letting a = 0.075 sample separate at 
T = 27.5°C. The resulting micellar volume fraction, calculated from the volume ratio of 
the phases and assuming a 2 % micelles in the lower N phase (see Fig. |2J) was = 0.14. 
Dynamic light scattering measurements on the upper isotropic phase are possible because 
the time for the two coexisting phases to remix is much larger than the time to perform 
the experiment. The micellar diffusion coefficient obtained in this way is compared in Fig 
[HI] with the Dm measured in sigle phase 5DHfiem. The overlap of the data confirms the 
absence of significant concentration dependence of D M at least within our experimental 
error. Moreover it supports the single component approximation, previously assumed in 
drawing the phase diagram (see Section Hi B|) . 

The deviation of Dm{T) from kBT/r)(T) indicates that the product kr,(f) , constant for 
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T > 40°C, grows appreciably as the system approaches the demixing transition. Altogether, 
the matching of Dm(T) with kBT/rj(T) at high T, and the independence of ko4> on , 
strongly suggests that, at high enough temperature, kr> ~ 0. Vanishing virial coefficients 
are not at all a surprise for water in oil microemulsions. Systematic investigations in systems 
prepared with different choices of oil and surfactants have revealed that inverted micelles 
may interact with a wide range of variable strengths, ranging from negative - hard sphere like 



- to positive virial coefficient 271] . This is reflected in other studies too. In one of them 28], 
it was found that a large body of data could be satisfactorily understood on the basis of hard 
sphere repulsion, with the hard sphere radius few A smaller than the measured hydrodynamic 



radius. Other investigations 
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301 ] concluded that contact attractive interactions had to 



be included to interpret the data. In another study 3_lJ it was actually found that collective 
diffusion was independent on droplet concentration, as in the 5DHfiem case here. The 
existence of attractive interactions in inverted micelles is usually attributed to an entropic 
gain coming into play when the hydrophobic tail of the surfactants of two droplets overlap 



32| . Our data indicate that such a " stickiness" interaction may be present in our system to 
compensate, as far as kn is concerned, for the hard sphere repulsion of the hydrophilic cores, 
certainly present. Thermal dependence of such interpenetration interaction also explains the 
inverted consolution curve ("lower" critical solution temperature) observed in the water- 



AOT-octane 
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]. We stress that such inverted demixing (induced upon increasing T) is 
the only gas-liquid like phase separation observed in water-in-oil microemulsions and that 
its reversed nature clearly demonstrates that the demixing observed in 5DHfiem (induced 
upon lowering T) cannot be explained by the same attractive potential. 



Previous works on DDAB 



focused on w/o microemulsions prepared by using hydro- 



carbon oils and with a much higher water/DDAB concentration ratio, propose a calculation 
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scheme for the micellar size which would yield, in our case R ~ 15.5 — 19. 5A, depending 
on how extended the DDAB chains are. Care should be taken, though, since the radius 
may depend on the nature of the oil. By comparing phase diagrams obtained with different 
kind of oils, two interesting elements appear: (i) DDAB is insoluble in alkanes if no water 
is added, indicating the existence of a minimum curvature of the DDAB layer, which also 



appears to be dependent on the alkane chain length I34L 135| | : (ii) DDAB is soluble in toluene 
with no added water but it does not form micelles 3J|. DDAB is soluble in the isotropic 
phase of 5CB with no added water. However, DLS measurements have revealed that, at 
variance with the toluene system, DDAB in the isotropic phase of 5CB form micelles. Their 
radius, measured in a sample with 15 % DDAB in 5CB by a different DLS set-up equipped 
by a more powerful laser source, is R ~ 18A. This number should be considered as very ap- 
proximate because no analysis has been done on possible intermicellar interactions. Overall, 
the solubility properties of DDAB in 5CB appears to be quite different from that in both 
alkanes and toluene. 

A simple estimate of the 5DHfiem inverted micellar structure is as follows. By assuming 



an intermediate (not fully extended) DDAB hydrocarbon chain length of 14A 24J, it follows 
that the water core radius R\y ~ 6 A. Since the DDAB/water volume ratio extends to 
the micellar structure, should the DDAB not compenetrate with 5CB molecules but rather 
form a compact spherical shell, we would have micelles of 13. 5A. This fact suggests a very 
strong interfingering of DDAB and 5CB molecules, in agreement with what implied by the 
solubility of DDAB in 5CB and by the large curvature of DDAB in 5CB. Overall, the DDAB 
volume in every droplet is small, and corresponds to about 13-14 DDAB molecules. 
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4- Micellar concentration fluctuations 

The growth of (Im) and of D(T)rj(T) /ksT as T DM is approached could, in principle, 
be explained in a variety of manners. It could be due to the growth of a 5CB layer of 
radially anchored molecules on the micelles surface, thickening as a consequence of the 
increasing nematic correlation length. Alternatively, it could follow from restructuring of 
the spontaneously assembled water-DDAB droplets. Finally, it could reflect the onset of 
attractive intermicellar forces. We will argument in favor of this last explanation. 

The development of a radial nematic layer around the micelles could certainly bring 
about a slowing down of their diffusion, but it certainly would not increase their scattering 
cross section for light. Indeed, given the size of the particles - well within the Rayleigh 



scattering regime |37[ , isotropically anchored molecules cannot be distinguished from a layer 
of isotropically distributed unanchored molecules. 

Restructuring of the self-aggregated nanodroplet could take place but only allowing for 
non-spherical shapes, since the ratio [DDAB]/ [water] has to be maintained in any aggregate. 
In order to account for a growth of about a factor 3 in the scattered intensity, aggregates 
should increment their volume of about 3 times. This would lead to the formation of rod-like 
micelles whose diffusion coefficient can be explicitly calculated and results to be at least 4 
times larger than expected for spheres of constant radius. Since the observed slowing down 
is only of a factor 2, we don't see room for pursuing this explanation. 

Interparticle attraction, besides slowing down their diffusion, also give rise to concentra- 
tion fluctuations, in turn enhancing the amount of scattered light. This effect is described 
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by the second virial coefficient kj in the virial expansion 



Im = Wl + M)" 1 (14) 

where I M ,o is the scattered intensity measured if interactions were not present, k^ — 
kj — kf, where both ko and kf are suitable integrals of the pair potential 26[. While 
kj describes concentration fluctuations, kf takes into account hydrodynamic interactions. 
Attractive interactions correspond to negative coefficients. For short ranged potentials, kf 
has typically the same sign of kj but is smaller in amplitude (about half or less). {Im) 
data in Fig. ^] are extracted from (Ivv) data on the basis of the ratio of slow and fast 
dynamics in the intensity correlation functions. Thus, as in the case of Iyv , the values of 
{Im) cannot be compared between different samples. However, the fact that hp ~ at high 
temperature, suggests that also kj is negligible in the same temperature range, as found in 
other systems having vanishing k^, 26|. The growth of {Im) as T decreases indicates that 
kj is temperature dependent, ranging from its "high T" value kf T = 0, to a value at the 
demixing temperature kj ~ —4.3 in the case of <fi = 0.15 and a value ki ~ —9.3 in the case 
of = 0.075 (from the data in Fig. Q2J). 

Since an increment in the attractive potential implies an increment in the negative value 
of both kj and kf, with the second smaller in modulus than the first, we expect kp to 
be negative and ko > kj. This prediction is in agreement with the experimental result, 
in which the effects of interactions are more conspicuous in the scattered intensity than 
in the diffusion coefficient. Therefore, while neither the growth of a nematic layer around 
the droplets nor micellar restructuring results compatible with experimental results, the 
appearance of intermicellar attraction is capable to account for the data. 

Although the specific nature of intermicellar interactions is not experimentally accessible, 
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their occurrence in the proximity of the demixing transition, in a temperature interval in 
which orientational fluctuations also develop, suggests the existence of a connection between 
the two phenomena. It is exactly to explore such a connection that we have developed the 
model described in the following Section. 



IV. THEORETICAL MODEL 



The nematic phase can be modelled in a variety of ways, having various degrees of realism. 
Striking, for its simplicity and effectiveness in capturing the strength of the first order I- 



N phase transition, is the Lebwhol-Lasher (LL) model 



a spin lattice model with the 



simplest orientational coupling compatible with the non-polar nematic symmetry. This 



model is also interesting because it enab 



a random field is properly added 
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es studying the effects of quenched disorder once 



40l | . Hence, the LL model appears as an attractive 



starting base to investigate the effect of nanoimpurities. The LL model can also be reckoned 



as the lattice counterpart of the widely known Maier-Saupe model 



repeately studied ever since 



431. 



42j, which has been 



The analysis of the pretransitional behavior of the pure LL indicates that, if a mapping 
is attempted between spins and molecules, a single spin represents not a single molecule but 
rather a cluster of about 50 molecules 



44|. In the case of 5DHfiem, the droplets in the 
solution are comparable in sizes with the length of the LC molecules, while filling the volume 
of 20-40 5CB molecules. Hence, one is then led to suspect that a lattice model containing 
a mixture of LC and micelles, taken to be of the same size, might be sufficient to capture 
the main qualitative features of the experiments presented so far. We now show that this is 



indeed the case. In our model, LC are represented by pointless unit vectors s which are free 
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to rotate in three dimensional space, but are fixed in position to a lattice of spacing a: 

s = (s , s 2 , s 3 ) = (sin 9 cos 0, sin # sin <f>, cos 0) (15) 

Micelles, on the other hand, are simply mimicked by vacancies, and are then represented 
- in the spirit of a lattice gas model - by a lattice variable a = 0(, 1) when the lattice site is 
vacant (occupied by LC). 

As usual, it proves convenient to introduce the traceless quadrupole tensor 

Q^(s) = 1(3^-0 a, (3= 1,2,3 (16) 

The final Hamiltonian then reads: 

-m = K qq e E q?q?°**v + K °° E E a ^y + E E h fQf** - E ^ ( 17 ) 

a/3 (xy) a/3 (xy) a/3 x x 

Here the first term reflects the fact that two LC have the general tendency to align 
along a common direction with coupling K qq = f3J qq , {(3 = 1/fcgT), the second term is a 
standard lattice gas Ising-like interaction with coupling K au = f3J aa , and the next two terms 
correspond to the external fields (h a ^ for the LC and fi for the micelles) controlling the LC 
orientation and micelle concentrations respectively. 

In the following the average fraction of micelles will be denoted as <fi. It is worthwhile 
noticing that in Eq (|T7|) the micelles interact both directly (through a K aa coupling) and 
indirectly, because of the presence of the LC. As we shall show, the latter are responsible 
for the most interesting features of the mixture, and are present even in the absence of the 
former. 
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A. Mean Field Analysis 



The diluted LL model has been studied using Monte Carlo simulations 45|, |46| In the 
following we shall study it within a mean field approximation akin that of its continuum 



counterpart 

To this aim we define a trial Hamiltonian 



-0H O = J2H bfQfa x + £ X x a x (18) 

a/3 x x 

where b al3 and A are variational parameters for the LC and micelles respectively. As 



remarked, at variance of Ref. 



8j, we shall exploit an Hamiltonian which is not simply 



factorized in two terms which depend on one variable only. Clearly a more general choice 



if = 5>({Qa,}) 



(19) 



(h ({Qcr x }) being a single site hamiltonian) could be exploited, but the calculations quickly 
become analytically unmanageable. 

Accordingly we shall use the following order parameters 



qf = (Qfv*) (20) 
n x = (cr x ) 

where averages are over the trial Hamiltonian. Note that <fi = 1 — n. 
The calculation now follows standard procedures. Once the mean field free energy /3Fmf 
has been obtained, it proves convenient to Legendre transform to a mixed function 

(3T MF = PF MF + hfqf + MX (21) 
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where /x* = /3/x. This function clearly depends upon the variational parameters b a ° and 
A, whose explicit values are obtained by the requirement that the free energy is stationary 
at those values. The final results for this mixed free energy in terms of the homogeneous 
value of the micelles order parameter n and of the LC order parameter q as identified by the 
diagonal matrix 



2 

2 












(22) 







q 



is 



(3r MF (q,n) 
N 



-K qq q 2 - 3K aa n 2 + -K qq qn + 



(23) 



+n logn + (1 — n) log (1 — n) — log47r — n log J ( — K qq q 




where we have defined the integrals 




(24) 



1. The phase diagram 



The isotropic-nematic transition is obtained by the conditions (when h a ^ = 0) 



d_ (3Y MF 
dq N 







(25) 



d_ f3T mf 
dn N 
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that is: 



( Y Kqqq 



(26) 



fj,* = 6K aa n - -K qq q - lo£ 



n 



1 — n 



log 



( Y Kqqq 



(27) 



where 



*(ah4^ /2(a) 



(28) 



2 2/ (A) 

Note that the Eq. (J27|) is a self-consistency equation for the nematic order parameter q. 

It is instructive to compare the above results with the corresponding one obtained in 
the absence of micelles, i.e. the fully occupied Lebwohl-Lasher model (FLL), and with the 
simple lattice gas case. The free energy for the former reads 



pt MF (q) 
N 



9 9 

-K qq q 2 + -K qq q- log(4vr) - log 



FLL 



[ Y Kqqq 



where the order parameter satisfies 



[q]fll = $ (^K qq q\ 



For the simple lattice gas (LG) case, on the other hand, we have 



(29) 



(30) 



(3t M F{n) 
N 



3K aa n(l — n) + n logn + (1 — n) log(l — n) 



LG 



(31) 



LG = U-^<7<T 



6K nrr n - 3K„ n - Ioe 



n 
1 — n 



(32) 
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A couple of comments are here in order. First the n — > 1 limit of Eq. (|22|) results into Eq. 
as it should. Note that since — 1/2 < Q aa < 1, the order parameter q is bounded and 
hence /i* — > oo yields n — > 1. Second Eq. (|24j) does not appear to be a simple combination 
of Eq. ()29|) and ()31|) as it has been argued on the basis of phenomenological theories 47 1. 

The final step necessary to obtain the coexistence values and nj (the nematic and 
isotropic LC densities at coexistence respectively) is the Maxwell equal area construction. 
It is clear that the resulting phase diagram will be influenced by the competition between 
the usual pure nematic-isotropic transition (associated to the free energy Eq. (|29j) ) and the 
usual lattice gas transition (whose free energy is Eq.(|3T)0. We shall study the isotropic- 
nematic transition as a function of the micelle molar fraction = 1 — n and on the ratio 
Jaa/Jqq- The results are depicted in Figs. EH EH] in the case J aa jJqq = (no direct 
micelle-micelle interaction), J aa /J qq = 3/4 < 1 (weak direct micelle- micelle interaction) and 
J aa j 1 J qq = 3/2 > 1 (strong direct micelle- micelle interaction) respectively. The temperature 
axis is scaled in terms of the nematic-isotropic (NI) transition temperature Tjvj of the fully 
occupied counterpart (FLL), which turns out to be 



Jqq 



1.9813... (33) 



In Fig. El there is no direct coupling among the micelles (J qq = 0) although there is an 
effective interaction as we shall see later on. As the micelles concentration increases, the 
isotropic phase (I) is depressed until it vanishes at = 1, a nematic re-entrant phase (N) 
appears in the small region, and the two phases are separated by a nematic-isotropic region 
(NI). In figure El the direct micelle- micelle coupling is weaker that the LC counterpart and 
the liquid-gas transition lies well inside the isotropic-nematic coexistence phase. Hence, the 
phase diagram is topologically similar to the previous one. At sufficiently high values of the 
direct coupling (Fig. ITol lower panel) the liquid-gas transition characteristic of the lattice gas 
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model preempts the isotropic-nematic transition at intermediate micelle molar fraction. In 
this case, the consolution curve separating an isotropic phase with low micellar concentration 
(IL) from an isotropic phase with high micellar concentration (IH), is symmetric because of 
the vacancy/filled site symmetry, and its highest point corresponds to the liquid-gas critical 
temperature T c 

k B T c 3 J (j (j 



/ 2 7 

°qq A 'Jqq 



The above phase diagram, as derived from a mean-field theory, is 



completely consistent 



471 ] . although the nu- 



with recent Monte Carlo numerical simulations on the same model 
merical values are not, as one expects from a mean- field-like theory such as the one treated 
here. 



2. Correlation Functions 

In order to make a connection with the experimental results probing both the paranematic 
fluctuations and the micellar density fluctuations, in this Section we compute the correlation 
lengths. This can be done rather conveniently within our mean field theory by computing 
the correlation functions as follows. The micelle-micelle correlation function is defined as 



97 y = ((1 " " *v)> " ((1 " *»)><(! - "»)> (35) 
= (o-xO-y) ~ {<Tx)(a v ) (36) 



This can be easily computed in momentum space 



c = ^E eip(x ^ CT<T (p) ( 3T ) 



N 

p 
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as its Fourier transform displays the following very simple form 



f° (p) = —77, ; ; n „ ^3 ( 38 ) 



n~ x {\ - n)" 1 - 2K aa J2a=i cos Pm 



a 



The correlation length for the micellar density £m can be immediately inferred from the 
behavior of g acr (p) in the continuum limit, leading to 



I e -\x-y\/^M 

9x°y = ~ A 7? 717 (39) 

with 



t_M_ I iW9 ( s 

a y n-^l-n)- 1 -6K aa 1 J 

It is worth noticing that £m — > as J^o- — > 0, meaning that, at the mean field level, 

there is no micelle-micelle interaction, in the absence of an explicit coupling among them. 

Accordingly, the above expression diverges at the (reduced) spinodal temperature Ts 

^ = 6n( l_ n )^ (41) 

Jqq Jqq 

which coincides with that of a simple lattice gas model. Only going beyond the mean field 
approximation, as described in the next Section, a LC mediated intermicellar interaction 
can be described. 

It is also interesting to evaluate how the micelle-LC coupling intrinsically present in this 
system affects the orientational pretransitional behavior. The nematic correlation function 
is defined as 

= (a x Qfa y Q^) - (a x Qf) (a y Q^) (42) 
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which can be conveniently expressed in terms of the derivative of the Legendre transform 
of F M f 



[G 



a/3fj,ul — l 



a.3 o au 



MF\ 



(43) 



dq x dqfy 

In the frame in which Eq.()22|) holds, it is then easy to calculate G 3333 = G xy . This can 
be computed along the same lines as for the micellar case. One obtains, after some algebra, 
an expression similar to Eq.(|39|) with the orientational correlation length £jv given by 



where 



3B(K, 



qq) 



A(K qq ,q,n)-3B(K qq ) 



(44) 



A(K qq ,q,n) = 12K qq {l-2K, 



L 99 



9 27 81 27 

— n(l - n)$ 2 (—K qg q) + — n R A {—K qq q) 



(45) 



B(K qq ) = 2K t 



</'/ 



(46) 



and where we have defined the quantity 



In particular in the isotropic phase 



h(X) /|(A) 
io(A) / 2 (A) 

: 0) this yields 



(47) 



6v 
a 



, i _ Ik v 
g =o y 5 !? 



diverging at the temperature 



A;rT* 9 



^qq 5 



(48) 



(49) 
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and reducing, in the limit n — > 1, to the correct value of the FLL. The line depicting 

the divergence temperature for the paranematic correlation length is shown in Fig. El as a 
function of the micellar density. 

B. High temperature expansion 

As previously remarked, one important question arising from the mean-field calculation 
above is whether there exists an effective micelle-micelle interaction, even in the absence of a 
direct coupling, which is induced by the fluctuation of LC degrees of freedom. One possible 
indicator of this is the presence of a micelle spinodal line which signals the approach of a 
demixing instability, i.e. the divergence of the correlation length, and we saw that there is 
no spinodal line in the absence on an explicit coupling among the micelles. We attribute 
this to a deficiency of the mean field calculation. 

To see that this is the case, we proceed as follows. First we split the original Hamiltonian 
Eq. lfTTj) at zero field (h a ^ = 0) in two parts 

H ({Q, a}) = H LC ({Q, a}) + H LG ({a}) (50) 

where 

- (3H LC ({Q, a}) = K qq £ E QfQFwv ( 51 ) 

a/3 (xy) 

and where the remaining term Hlg {W}) is a lattice gas hamiltonian independent on the 
LC degrees of freedom {Q(s)}. We then introduce an effective micelle-micelle reduced 
Hamiltonian — j3H e ff({a}) by integrating out the LC degrees of freedom of the original 
Hamiltonian EqloTH 
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where we have defined 



Ds 



n^=n 



2jx 



+ 1 



d(cos9 x 



(53) 



In the high temperature limit T » T c then K qq << 1 and thus 



Ds e Kqq ^ a @^( xy '>®* l3 ®y a ' Tx ' Jv 



(54) 



Ds 



i + K qq e E Qfotrwv + \ K l E E Qfotrwv E E oro2r^« + • • • 



a/3 (13/) 

Using the following results 



a/3 (xy> 



[s) = 



(55) 



J dsQ aP (s)Q» v (s) = -\<k^^ v + \tz{8°^S Pv + cP^) 



(56) 



5 5 

we obtain, to first order in K qq and after rexponentiating the result, a new lattice gas 
Hamiltonian 



-(3Hef f (M) = K e JJ E + E M ~ a ^ 



where 
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K eff = K H K 2 

which amounts to an effective coupling 



9 J 2 

T e ff = T 4- qq 
J ™ + 40 k B T 



(57) 



(58) 



(59) 
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Notice that J^J depends upon the temperature T. 

In the case K au = 0, we can calculate the spinodal line Ts(n) on the basis of the mean 
field solution of the standard lattice gas model. We obtain 



k B T s 27 , 

if-vs"* 1 -) (60 » 

which is also depicted in Fig. Note that this result differes from that of a simple 
lattice gas (Eq. I3T|) because of the T dependence of the effective coupling. 

This analysis shows effective micelle-micelle interactions, can be computed by eliminating 
the LC degrees of freedom from the problem, as done here through a calculation correct to 
a first order in a high temperature expansion. 

C. Radial coupling 

Before ending this section, we wish to comment on a possible additional term appearing 
in the general Hamiltonial (|17|). As described in Ref. 7(, the properties of the present system 
can be strongly influenced by a directional coupling between the LC and the micelles. Hence 
one is tempted to include a tendency of the LC to be oriented along the radial direction 
of the micelles whenever the LC lie in their proximity, as a possible consequence of the 
tendency for homeotropic anchoring of DDAB. 

To account for this we add the following additional term —/3AHa to the Hamiltonian in 
Eq.fUZD 



-P*Ha = frEE^-^K^f (y-x) Q (y-xf +a x (l-^)Qf (x-yHx-yf] (61) 

a/3 (xy) 

As described in Fig. the sign of K qa = /3J qcr identifies the preferred orientation of the 
LC with respect to its distance from the center of the neighboring the micelle. If J qa > 
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the LC at the LC-micelles interface tend to align radially. It is easy to see that, at the mean 
field level, the above term does not affect the phase diagram. This is because mean field 
theory effectively eliminates the site dependence of the order parameters. 

Things are different when the same coupling is investigated through high temperature 
expansion. In this case, the radial coupling leads to non-trivial three-body terms. This, at 
variance with the mean-field theory, indicates the possible relevance of the LC anchoring on 
the mi cellar surface. 



V. DISCUSSION 

A. Comparison between experimental results and theoretical model 

The existence of Casimir-style, fluctuation mediated inter-micellar attractions, demon- 
strated by the high temperature expansion in the case of no direct coupling [K aa = 0), is 
a key result that enables us to interpret the observation of pretransitional micellar density 
fluctuations in systems where paranematic fluctuations are also observed. Indeed, the theo- 
retical modelling enables us to understand the inter-micellar attraction documented in the 
experiments as interactions mediated by paranematic fluctuations. 

More specifically, we argue that the experimental conditions we investigated correspond 
to K aa w and K qa rs 0. The simultaneous presence of micellar and orientational fluc- 
tuations, and the fact that upper critical solution points have never been observed in w/o 
microemulsions, make us reject as unlikely (but not impossible) that a direct interaction 
may be at the origin of the observed behavior. Such an interaction should in fact be in- 
dependent from the ordering of the solvent and of a different nature with respect to those 



leading to lower critical consolution point in other systems 



3311 . On the other hand, the 
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fact that the single molecule dynamics detected by EBS is unaffected by the presence of the 
micelles (see Section |IHB2)) . together with the nanosize of the droplet radius, way shorter 
than typical extrapolation lengths for surface anchoring, strongly support the notion that 
the 5CB anchoring at the micellar surface is very weak. 

Therefore, we read the observed increase in (I M ) and decrease in D M /D on the basis 
of the theoretical predictions that the softening of the micellar system (increase of micellar 
osmotic compressibility and associated slowing down) is expected on the only basis of the 
micellar dilution effect of the local nematic matrix. The existence of a bell-shaped spinodal 
curve in Fig. O suggests that a similar curve could be present in the experimental phase 
diagram as well (Fig. |2J, although hidden into the coexistence region. Thus, on lowering T, 
i.e. on approaching the spinodal line increasing micellar density fluctuation are expected. 
Accordingly, we have fitted (Im) vs. T with (T — Tg) _1 . The extracted values for Ts(4>), 
which suffer of a rather large indetermination because of the intrinsically limited range of 
the fit, are shown in Fig. |21 The two data point appear to decrease with increasing for 
(f) > 0.075, but for small experimentally inaccessible <p, Ts(4>) must take a positive slope 
since the system of highly diluted micelles necessarily approaches gas-like non-interacting 
behavior. 

The comparison between the experimental phase diagram (Fig. |2J) and the theoretical 
one (FigEJ) enlightens a remarkable similarity in all the qualitative features: re-entrant 
nematic phase, shape of the coexistence curve, non-trivial divergence lines determining the 
orientational and micellar pretransitional fluctuations. However, no quantitative comparison 
can be drawn between the two phase diagrams because of the approximations at the base 
of the model and in the techniques used to solve it. Both the temperature and the micellar 
density axes cannot be easily mapped on the corresponding experimental quantities. For 
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example, the T^i — T* gap for the pure nematic as well as the value of corresponding to 
the maximum in the spinodal curve are both much larger than the corresponding observed 
quantities. These shortcomings are consequences, respectively, of the intrinsic limitations of 
the mean field approximation and of having assumed that micelles and LC molecules have 
equal volumes . 

According to the interpretation here proposed, the expulsion of micelles by orientational 
ordering, explicitly manifested at low T in the macroscopic phase separation, is at the root 
of the whole phase behavior of 5DHfiem. Nematic correlations appear more easily where 
the local micelle concentration is lower, and in turn force out micelles, thus further reducing 
the local amount of annealed disorder and increasing the local T*, with the consequent 
increase of local nematic order. Theoretical considerations indicate that sufficiently small 
impurities incorporated in the nematic phase display attractive long ranged interactions 



due to orientational fluctuations of the LC 



481 ] . This is a consequence of the impurity 



induced depression of the local nematic order which are thus segregated to minimize the 



total free energy, a mechanism that bears some resemblance to the 



colloidal aggregation 



jrewetting induced 
We argue that 



or to the depletion forces in critical systems 
the same mechanism should be present in the isotropic phase but with an exponential cut- 
off at £at. This notion is corroborated by the experimental coincidence between the average 
intermicellar distance d and the maximum value of £,N,max = £o[T*/(Tdm — T*)]^, with 
~ 10A : for = 0.075, £ NtTnax ~ 95A, d ~ 58A; for = 0.15, £ N , max ~ 75A, d ~ 74 A). 



B. Concentration dependence of the transition temperatures 

All the observed 5DHfiem behavior appears determined by three relevant temperatures: 
the demixing temperature T DM ((/)), the divergence temperature for paranematic fluctua- 
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tions T*(<j)), the spinodal temperature Tg((j)). We want here discuss the origin of their 
dependence. 

One of the most striking features of the phase diagram of the 5DHfiem is the high 
negative slope of T DM (</>) and T* ((/)): T NI - T DM {<f)) ~ 4.4 A and T NI - T*(</>) ~ 7 A for 
the 5M/iem <p = 7.5%; T NI - T DM (<j>) ~ 7.6K and T NI - T*{4>) ~ 11 A' for the 5DHfxem 
4> = 15%. Thus dTr>M(<P) / d(p ~ — 60A" and dT* ((j>) / dcj) ~ — 80A'. These numbers are quite 
larger than the T shift in the position of the specific heat (c p ) peak observed by studying 
LC with nanosized quenched disorder: dT NI ((fi)/d(j) ~ —7.1 A for 8CB embedded in silica 



aerogel 



and dT NI ((j))/d(f) ~ — 8.2A for 8CB incorporating soft gels of aerosil particles 



. The origin of such a difference will also be discussed in the following. 



1. Demixing temperature (Tjjm) 

The determination of the presence and size of the droplets enables us to better discuss the 
low concentration part of the phase diagram. Indeed, the value of j3i, too low if calculated on 
the basis of the molar fraction of monomeric DDAB (see Section III. A), takes a much higher 
value when referred to the molar fraction of micellar aggregates, which can be calculated 
from the mass of a single micelle. The radius that micelles would have if DDAB were 
compacted in a sphere was estimated to be about 13. 5A (see Sec. III.C.3). Based on this 
value, the slopes of the isotropic and nematic phase boundary vs. contaminant fraction 
are, respectively (3i ~ —4.17 and (3n ~ —31.85. We thus obtain (fi^ 1 — (3j l ) ~ 0.21, in a 
much better agreement with the prediction of ((3^ 1 — fij 1 ) ~ 0.25. The numerical difference 
between the two numbers indicates that the calculation partially underestimates the slope 
of the isotropic phase boundary. 

The same relationship between phase boundary and contaminant concentration can be 
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reformulated under a different prospective, which offers a better insight on its physical basis. 
By approximating (3jj l ~ 0, the occurrence of the 5DH[xem demixing transition coincides 
with a transformation of a fraction of the LC from isotropic to nematic and, at the same 
time, with the segregation of the inverted micelles in the isotropic phase. By balancing the 
gain in free energy density for the nematic ordering of 5CB (AGni{T)) against the work 
done against the osmotic pressure II of the micelles, we obtain 



AGW(T) = 11(0) (62) 
AGni{T) can be computed on the basis of the Landau-deGennes expansion, whose coeffi- 



cients for bulk 5CB have been determined 



Indeed, since the micelle-disordered isotropic 



cannot differ much from the unperturbed isotropic (especially because surface coupling is 
very weak), and since the nematic phase in the 5DHfiem is almost micelle-free, AGjyi{T) 
calculated for bulk 5CB is a good approximation for the difference in free energy between 
nematic micelle poor phase and isotropic micelle rich one. The osmotic pressure, on the 
other side, can be calculated on the basis of and R, the micellar radius, and by assuming a 
perfect gas behavior. With this with these ingredients, it is possible to extract fli in the limit 
of small . With this crude estimate we obtain (3j ~ —6.5. Although overestimating the 
measured j3i, this figure is close enough to the data to indicate that indeed the temperature 
depression of the demixing transition can be understood as a direct consequence of the very 
nature of the demixing. This fact once again confirms the appropriateness of describing the 
hDHfiem as a pseudo-single component system. 



14 



2. Paranematic divergence temperature (T*) 

According to the model T*(0) = T*(0 = 0) ■ (1 — 0) (see Section HVJ) . The decrease of 
T* with can be easily interpreted as a dilution effect on the coupling strength. Indeed, 
the dilution brings about a decrement in the average coordination number of each spin. As 
a consequence J qq , the coupling coefficient averaged on all bonds, also decreases. Since is 
the probability to find a micelles as nearest neighbor of each given spin, J qq — J qq ■ (1 — 0). 
Hence, it appears that the easiest interpretation of the observed T*(0) is to picture it 
as a direct consequence of a dilution induced weakening of the molecular field, which, in 
the experimental case, is associated to the presence of interfaces reducing the coordination 
number of the 5CB molecules. A possible estimate of the fraction of 5CB molecules at contact 
with the micellar surface is given by the volume between 13. 5A and 20A from the center of the 
micelle, i.e. by the difference between the volume hydro dynamically attached to the micelle 
and the total water-DDAB volume. This would imply that a fraction of about / = 3.2 • 
of the total 5CB molecules experience reduced coordination because of the proximity of 
DDAB molecules. Comparing T*(0) = T*(0 = 0) x (1 — 0) with the experimentally observed 
dT* (0) / d(j) ~ — 80K, suggests that the coordination number of 5CB molecules at the DDAB 
interface is reduced of 8 %. This number appears reasonable, since we expect it to be lower 
than the simplest available analogy of a hard sphere crystal limited by a plane. In that 
case spheres facing the plane have a coordination of 9 while in the bulk the coordination is 
12, yielding a decrease of 25% of the coordination number. This simple evaluation, despite 
its extreme roughness, supports the notion that the concentration dependence of T*, and 
thus the whole T downshift in the pretransitional paranematic behavior, is due to dilution 
induced weakening of the intermolecular coupling. 
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3. Comparison with LC in silica gels 

The N-I transition in LC with dispersed silica nanostructures, either aerogel or aerosil, 
displays a specific heat c p peak downshifted in T with respect to the corresponding transition 
in the bulk and downshifted in T. In those systems it is found that the whole paranematic 
behavior is unmodified but "rigidely" shifted in temperature 



). It is quite remark- 
able that such observed dT* /dip in the aerogel/ aerosil LC systems is about ten times smaller 
than in 5DHfiem. The different nature of the disordering provided by the two structures 
hosted in the LC, mobile "annealed" disorder in 5DHfiem and static "quenched" disorder 
in nanosilica-LC materials, cannot be the origin of this large effect. This is because the dif- 
ferent degrees of freedom are highly mobile in the high temperature phase, thus preventing 
a clear distinction between a frozen impurity and an equilibrated one. 

The difference between the two systems can, at least in part, be understood on the basis 
of different dilution induced weakening of intermolecular coupling, along the line of the 
discussion in the previous section. X-ray characterization of the empty aerogel structure - 
not available in the aerosil case - enable extracting the specific interface extension of the 
structures, expressed in the S/V ratio, where S is the silica surface area and V the aerogel 



For a particular family of silica aerogel extensively characterized in the past, it is 

i 



volume 

found that, approximately, S/V = 0.1 lA 



It follows that the fraction of LC molecules 
at the silica interface is approximately / ~ S/V ■ £/vlc ~ 0.8 , where Vlc is the volume 
of a single LC molecule. Thus, at equal volume fraction, the silica aerogel is characterized 
by an interface about three times smaller than the 5DHnem, which would imply a value of 
dT*/d(j) three times smaller, thus only in part accounting for the experimentally observed 
difference. This fact indicates that some other factor plays a relevant role in determining 
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the T* shift, possibly related either to the different surface coupling of the two surfaces or 
to the different geometry of the disordering hosts: well dispersed in the 5DHfiem, while 
necessarily interconnected in the aerogel. 

VI. CONCLUSIONS 

In this article we have presented a systematic investigation of a system where thermody- 
namically stable nanoparticles are dispersed in a thermotropic liquid crystal. As it appears 
from the experiments, the self aggregated nature of the particles does not play any observ- 
able role, in agreement with the notion that such self assembly is governed by forces larger 
than the molecular field strength provided by nematic ordering. This system enables study- 
ing the effect of annealed disorder on the phase diagram of the liquid crystal, and, at the 
same time, thanks to the large optical contrast provided by the droplets' water core, it also 
allows direct observation of the statistical properties of the dispersed nanoparticles. 

We found that the system exhibits a phase transition combining the demixing of the 
nano-droplets with nematic ordering in one of the coexisting phases. The phase diagram 
also features an intriguing reentrant nematic phase. This feature indicates a temperature 
dependent solubility of nano-droplets in the nematic phase. 

The characteristic frequency of the electric birefringence spectra expressing the time to 
redirect the molecular axis, shows that the dispersed nanoparticles do not alter appreciably 
the local molecular environment which remains bulk-like. This fact supports the notion that, 
mostly because of their nano-size, microemulsion droplets do not provide significant anchor- 
ing on the neighboring LC molecules. Hence, the observed effects of the nanoinpurities on 
the IN phase transition have to be understood simply as consequences of dilution, in strong 
analogy with the most conventional understanding of the effects of molecular contaminants. 
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Along this line we have here presented a theoretical analysis aimed to determine the phase 
behavior of a diluted Lebwohl-Lasher spin lattice model, where LC molecules are mimicked 
by pointless spin vectors on a cubic lattice and micelles are vacancies on the same lattice. 
The spin variables interact via a Maier-Saupe orientational coupling, provided that nearest 
neighbour pairs are occupied by LC variables. This system has been solved by mean-field 
theory and high-temperature expansion. 

Experimental observations and theoretical predictions nicely match offering a simple 
ground to interpret all the basic features of the system. 

The rather complex topology of the observed phase diagram (demixing transition and 
re-entrant nematic phase) simply follows from the tendency of liquid crystal towards ori- 
entational order with no need for either micelle-micelle or micelle-LC direct interactions. 
The dependence of phase boundaries and divergence temperatures on the concentration of 
impurities conveys the mutual effects of orientational order and nanoparticle dispersion. 
Furthermore, the study of the light scattered by the nanodroplets in the isotropic phase 
close to the IN transition, enlightens the presence of intermicellar attractive forces which 
the theoretical model enables interpreting as a consequence of the paranematic orientational 
fluctuations. This can be seen by eliminating the LC variables in the Hamiltonian of the spin 
model, approximated to the leading order in a high temperature expansion. An effective 
temperature dependent micelle-micelle interaction is found, whose strength depends on the 
orientational coupling between the LC molecules. This fluctuation mediated interaction is 
intrinsic to the very presence of nanoimpurities, necessarily diluting the LC, while the ad- 
dition of possible surface anchoring does not change the picture within a mean-field theory, 
but has a non-trivial effect within our high-temperature expansion. Accordingly, it appears 
interesting to explore how the phase diagram is modified by increasing the droplet size 
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and the hydrophobic portion of the surfactants, which should lead to progressively stronger 
anchoring. 

We acknowledge support from NATO (Collaborative Linkage Grant) and from MIUR 
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Sample 4"=0 DS 




00 (b| 

FIG. 1: Pictures showing the phase separation of 5DHfiem. At low temperature (a) the system 
consists of a nematic phase filling 31 % of the total volume, coexisting with a lighter isotropic 
phase. At higher temperature (b) the system appears as a single isotropic phase. 
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0.0 0.1 0.2 0.3 0.4 0.5 

FIG. 2: Experimental phase diagram featuring isotropic (I), nematic (N, grey shadowing) and 
phase coexistence (I+N). Full dots: phase boundaries obtained from the volumetric measurements 
shown in the inset. Open dots: Tjjm of the four samples used for volumetric measurements. Open 
diamond: extrapolated divergence temperature T*(</>) of paranematic fluctuations. Open squares: 
extrapolated divergence temperature Ts{4>) for micellar density fluctuations. Inset: volume fraction 
of the isotropic phase in the two phase coexistence region for: (f> = 0.02 (squares), (f> = 0.05 (circles), 
(f> = 0.075 (diamonds) and = 0.15 (triangles). 
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FIG. 3: Electric birefringence spectra measured on a sample of bulk 5CB at various temperatures: 
T = 35.93°C (triangles), T = 35.53°C (dots), T = 35.28°C (squares), T = 35.11°C (diamonds). 
Lines: fittings with the Debye Eq. |1J Inset: enlargement of the high frequency data and fits. 
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FIG. 4: Electric birefringence spectra measured on a sample of 5DHfj,em with <f> = 0.15 at various 
temperatures: T = 33.66°C (triangles), T = 32.86°C (dots), T = 31.85°C (squares), T = 30.80°C 
(diamonds). Lines: fittings with the Debye Eq. 0J 
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T[°C] 

FIG. 5: Scattered intensity in in a 5DHfiem sample with <f> = 0.15 (main figure) and in bulk 5CB 
(inset) as a function of temperature. Empty dots: depolarized scattered intensity {Ivh)- Full dots: 
polarized scattered intensity (iyy). Crosses: depolarized scattered intensity scaled by a 4/3 factor 
(see Eq. |SJ). Lines: fittings with (T - T*) _1 . Main figure: T* = 21.9°C, indicated by the vertical 
line. Inset: T* = 33.8°C, for both the polarized and depolarized components. 
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FIG. 6: Comparison of the depolarized intensity scattered by bulk 5CB (circles) and bDHfiem 
with (f> = 0.075 (squares) and (j) = 0.15 (diamonds) as a function of temperature. Lines: fittings 
with (T - T*)" 1 yielding T*{bulk) = 33.8°C (solid vertical line), T*((j> = 0.075) = 27.2°C (dashed 
vertical line) and T*(4> = 0.15) = 21.9°C (dotted vertical line). 
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FIG. 7: Normalized autocorrelation functions of the polarized (VV, full dots) and depolarized 
(VH, open dots) scattered intensities for the <p = 0.075 bDH^iem sample at T = 31°C. Inset: 
comparison of normalized autocorrelation functions of the polarized scattered intensity for the 
4> = 0.075 (full dots) and <6 = 0.15 (diamonds) hDH^em samples at T = 31° C. Lines: fit with 
double exponential decay (VV) and single exponential decay (VH) 
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FIG. 8: Low frequency Kerr constant B of Bulk 5CB (circles) and 5DHfj,em with <p = 0.075 
(squares) and 4> = 0.15 (diamonds) as a function of temperature. Crosses: high frequency Kerr 
constant asymptote B e of bulk 5CB multiplied by a factor 15 to overlap B&- Lines: fittings with 
(T-T*)- 1 , as detailed in the text, yielding T*{bulk) = 33.4°C (solid vertical line), T*{(j) = 0.075) = 



27.3°C (dashed vertical line) and T*(<j> = 0.15) = 21.8°C (dotted vertical line). 
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FIG. 9: Arrhenius plot of v c vs. 1000/T for bulk 5CB (circles) and for 5DH^iem sample with 
4> = 0.15 (diamonds). Dashed line: Arrhenius fit of all data. Vertical lines indicate Tjv/ (solid line) 
and Tdm{4> = 0.15) (dotted line). 
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FIG. 10: Comparison between characteristic times and inverse frequencies characterizing the dy- 
namic light scattering and electric birefringence spectra. Full dots: Tp^yv for bulk 5CB. Empty 
dots: tpn,vh for bulk 5CB. Squares: tpn,vh for the <j) = 0.075 5DH^iem. Diamonds: tpn,vh 
for the <p = 0.15 bDHfiem. Pluses: 1/(2ttz/ c ) for bulk 5CB. Crosses: l/{2%v c ) for the </> = 0.15 
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FIG. 11: Paranematic correlation time tpn divided by the shear viscosity rj for bulk 5CB (circles) 
and f>DHjiem with (ft = 0.075 (squares) and (ft = 0.15 (diamonds) as a function of temperature. 
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FIG. 12: Intensity scattered by microemulsion droplets {Im) extracted from dynamic light scat- 
tering data as a function of temperature. Squares: <ft = 0.075. Diamonds: eft = 0.15. Crosses: {Im) 
obtained from the difference between Iyv{4> = 0.15) and 4/3 • Iyn{§ = 0.15) and normalized to 
match diamonds at T = 60°C. Lines: fittings with (T - T s )~ l yielding T s (4> = 0.075) = 19.4°C 
(dashed vertical line), T s (4> = 0.15) = 14.6°C (dotted vertical line) 
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FIG. 13: Micellar diffusion coefficients Dm as a function of temperature for bDHfiem samples 
with (f> = 0.15 (full diamond), (f> = 0.075 (squares). Empty diamonds Dm measured in the upper 
phase of a = 0.075 5DHfiem at T = 27°C where independent estimates indicate a micellar 



concentration of <j) = 0.14. Line: expected temperature dependence of Dm for non interacting 



microemulsion droplets having radius R = 18A calculated on the bases of rj(T) from Ref. 
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FIG. 14: Theoretical phase diagram calculated in the case of no direct intermicellar interaction 
J aa = 0. The temperature axis is scaled with the nematic-isotropic transition temperature Tjvj 
of the fully occupied model (see Ea. Q33[) ). The dotted line is the spinodal line T$ for micellar 
demixing, obtained through a high temperature expansion, Ea. Q6UJI . The dashed line corresponds 
to the divergence temperature T* for the nematic correlation length, Ea.(|49|). 
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FIG. 15: Theoretical phase diagram in the case of direct intermicellar attraction of strength 
Jaa/Jqq = 3/4 (upper panel) and J a a/Jqq = 3/2 (lower panel), where the value corresponding 
to (j) = 0.5 is the demixing temperature Tq given by Eg. ((511 ) 
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FIG. 16: Possible Liquid crystal couplings at the micelle surface. Left: radial coupling. Right: 
tangential coupling, a: 2D pictorial view of the couplings, b: schematization of the coupling as in 
the lattice spin model, where circle represents an empty site and the segments the headless spins. 



coefficient K qa of eq. 1611 expressing the coupling strenght. 
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